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Abstract 



We consider density fluctuations of a two-fluid model consisting of hydrogen 
plasma and radiation prior to the cosmic hydrogen recombination. As investi- 
gation method that of the dispersion relations is applied, which have been de- 
rived from the general-relativistic sound-wave equations taking into account 
the coupling between plasma and radiation carefully. We obtain growing un- 
stable acoustic modes within the mass range 2 • 1O 6 M < M < 6 • 1O 12 M . 
In a second step the coupled differential equations for the amplitudes of the 
unstable modes are integrated numerically with respect to time where the 
integration extends from the initial time prior to the hydrogen recombination 
up to the present time. We find a significant increase of the amplitudes up to 
4 orders of magnitude, if the Universe is described by a cosmological model 
with a positive cosmological constant (A ~ 2, 2 ■ 1CT 56 cm -2 ) and a posi- 
tive curvature (Lemaitre-Universe) without an essential amount of cold dark 
matter. We conclude that the existence of galaxies confirm these statements. 
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1 Introduction 



In a previous paper (GeBner & Dehnen 2000) we have investigated stability 
and instability of the density fluctuations of a cosmic two-component fluid 
consisting of hydrogen plasma and radiation for explaining the galaxy forma- 
tion using Newtonian cosmology and the method of the dispersion relation. 
Although the investigation of cosmic density fluctuations has been performed 
already by several authors (c.f. the book of Coles & Lucchin 1995), the es- 
sential new idea in our paper was a careful treatment of the interaction of the 
plasma and the radiation prior to the hydrogen recombination taking into 
account that the radiation pressure of the photons on the plasma acts only 
on scales comparable or larger than the free path length of the photons (c.f. 
Diaz- Rivera & Dehnen 1998; Rose 1993, 1996). In this way we were able to 
explain the mass-spectrum of the galaxies in detail by the spectrum of the 
unstable fluctuation modes with a lower and an upper mass limit of about 
1O 6 M and 1O 13 M respectively. 

An open question is however the time evolution of the unstable fluctua- 
tions, especially its dependence on cosmological parameters as the cosmologi- 
cal constant or quintessence, the space-like curvature and dark matter. Until 
now it is commonly assumed, that most of the matter in the Universe exists 
in a nonbaryonic dark form, in order to accelerate structure formation on 
large scales (Peebles 1993, 2000; Hogan & Dalcanton 2000; Hu et al. 2000). 
But we will show, that such an assumption is not required. Furthermore it 
is interesting to see, how the Newtonian calculation of our previous paper 
follows from a general-relativistic one. 

Therefore in the present paper we formulate at first the fundamental 
equations general-relativistic and subsequently perform a first order pertur- 
bation calculation in order to derive the sound-wave equations for the plasma- 
(matter-) and radiation-densities. From these we proceed to the Newtonian 
limit and demonstrate that this step is justified so long as the masses of the 
matter fluctuations are smaller than 10 19 M Q (see also Bonnor 1957). For 
lower masses we determine using the method of the dispersion relation once 
more the unstable fluctuation modes and confirm the Newtonian results of 
our previous paper by the general-relativistic approach. 

In a second step we investigate the time evolution of the unstable fluctu- 
ations by a numerical integration of the sound-wave equations started from 
a point in time prior to the recombination of the hydrogen until the present 
time; the initial values for the numerical calculations are chosen as identical 
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with those of the unstable modes according to the dispersion relation. We 
obtain the following interesting result: The increase of the unstable fluctua- 
tion amplitudes is, as it is well known (see eg. Weinberg 1972), not described 
by an exponential function of time but rather by a power law, since the co- 
efficients of the sound-wave equations are time-dependent due to the cosmic 
expansion, and it is argued, that this increase is not sufficient for galaxy 
formation and that therefore the assumption of dark matter is unavoidable. 
However this situation changes drastically in favour of an exponential law 
of time with an increase of the fluctuation amplitudes up to 4 orders of 
magnitude, if a positive cosmological constant or quintessence is present and 
additionally the space-like curvature is positive! Therefore we conclude that 
the existence of galaxies confirm these last two facts in agreement with the 
results of W. Priester et al. (1995) and the new observation of an acceleration 
of the cosmic expansion at the present time (Riess et al. 1998). 

Furthermore, the mode structure of our two-fluid model depends essen- 
tially on the ratio of radiation and matter in the Universe. If this ratio 
differs only slightly from the used value, the modes and their time evolution 
change, which means that the observed structures in the Universe reflect im- 
mediately the ratio of radiation and matter. In this context it may be of 
interest, that additional cold dark matter does not improve the comparison 
with the observations in contradiction to the common opinion. 

2 The fundamental general-relativistic equa- 
tions for the interacting two-fluid model 

Prior to the recombination of the hydrogen we consider the electron-proton 
plasma and the photon gas as two components of the cosmic fluid. The undis- 
turbed homogeneous fluid is assumed to be in thermodynamic equilibrium. 
The matter density of the plasma is p p and that of the radiation (photon 
gas) p r = <tT 4 /c 2 (a Stefan-Boltzmann constant, T radiation temperature). 
The plasma 4- velocity is denoted by t> M (f M f M = 1) and that of the radiation 
flow by c M (c M c M = 1). Then the energy-momentum tensor for plasma and 
radiation is given by: 

(p) 

Tp, = (ppC 2 +P P + P^e + Pe^ ) V" ~ (pp + Py^c + Pc^-y ) 9iw ( 2 -l) 
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and 



Tfn, = (p r C 2 +p r + Pe^ 7 + P^_ e ) C V C " ~ (Pr + Pe-vy + P^-e) ( 2 ' 2 ) 

with the energy-momentum laws£] 



(P) (p«-r) ( r ) (r«-p) 

lit- = > hi/ = C/u > (2-3) 



where 



(p<-r) (rt-p) 4 c 2 y gC 

^ =-C^ = - pv -Tic,, - Vf,) + p p — {Cf, - Vf,) (2.4) 



O ("V 6fj 



describes the energy-momentum transfer between radiation and plasma (ki- 
netic energy head). Herein 



m 3 m 



^7 = _ j '•e = ~ _ \ 4-0) 



are the free path lengths of the photons (Z 7 ) with respect to the plasma and 
of the plasma particles (l e ) with respect to the photon gas. Insertion of (|2.5|) 
into (|2.4| ) shows that the second term of the right hand side of ( |2.4j ) can be 
neglected as long as the mean thermal velocity v e of the electrons is significant 
smaller than the velocity of light, i.e. v e <C c. This is in fact the case up 
to redshifts of the order of 10 8 , that means for our purpose. Of course, the 
interaction of the photons with the plasma takes place essentially via the 
electrons and therefore in ( |2.5| ) a Th = ^-(^f^j) is the Thomson scattering 
cross-section of the electrons. On the other hand, the inertia of the plasma 
is determined by the protons strongly connected with the electrons, so that 
in (|2.5|) m denotes the proton mass. 
Furthermore 

Pr = \PrC\ P P = P P C? With C s 2 = 2^ (2.6) 

represent the pressures of radiation and matter in (|2.1|) and (|2.2|) respectively 
and c s is the isothermal sound-velocity of the plasmaQ, which is considered 

1 H z/ means the covariant partial derivative with respect to the coordinate x u , and \ v 
the usual one. 

2 We prefer the isothermal sound- velocity because the plasma is embedded in the enor- 
mous heat-bath of the radiation. The ratio of the thermal energy density of radiation and 
matter amounts to 3 • 10 9 . 
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as an ideal Boltzmann gas. The remaining pressure-terms in ( |2.1| ) and 
describe the pressure interaction between plasma and radiation limited by 
the free path lengths (|2.5|) as follows: 



(a) p^ c = -p^ e = el L > ± p r , (b) p c ^ = -p e % = 4 p p , (2.7) 

where A signifies the covariant 3-dimensional Laplace-operator with respect 
to the rest frame of plasma and radiation. Equation (|2.7| a) represents the 
hydrostatic pressure of the photon gas on the plasma (radiation pressure) and 
its backreaction (R) on the photons; relation fl2.7|b) describes analogously the 
hydrostatic pressure of the plasma on the photon gas and its backreaction (R) 
on the plasma. A derivation of the expressions (|2.7|) is given in Appendix A. It 
is remarkable that the limitation factors in (|2.7| ) are operator- valued, so that 
an assumption of the special shape of the matter and radiation fluctuations 
is not necessary in the beginning. 

Finally, for the determination of the metric of space-time we use Einstein's 
field equations of gravitation with the cosmological constant A (G Newtonian 
gravitational constant): 

Rv» = -^T (T^-lTg^)+Ag^, (2.8) 



where according to ( 2.1|) , (2.2) and 



(p) CO 

T^ u = T^ u + T^v = p p c v^Vy - p p (g^ - v^v u ) + p T c c^c u - (2.9) 



with 

W = (2.10) 



guaranteed by ( |2.ci| ) and ( |2.4| ). Due to the last term of the right hand side of 
( |2.9| ) the reaction principle in a relativistic theory is not fulfilled exactly. 



3 Perturbation theory 



For the investigation of the cosmic matter and radiation fluctuations we start 
from the undisturbed Robert son- Walker Universes with the line-elements in 



6 



isotropic space- like coordinates 

ds 2 = dx 42 — — — ^5 ik dx l dx k , r 2 = 5 ik x t x k , x A = ct, (3.1) 

(1 + \er 2 ) 

where e = +1, 0, —1 is valid for the case of positive, vanishing and negative 
curvature of the 3- dimensional space t = const.. Subsequently we perform a 
first order perturbation ansatz as follows 

9i*> = 9»v + 9p,(x a ), v ft = Vn+Vn(x a ), c M = c M + c M (x a ), 

(3-2) 

Pp =Pp(t) +Pp(B a ), Pr =Pr{t) +Pr(x a ) 

with 

>=>=$»= c„ = (0,0,0,l) (3.3) 
and with respect to the normalization of t> M and c^: 

1 4 1 4 1 1 1 1 to A\ 

v = c = — v± = -c 4 = — g 9m- (3.4) 

Herein we consider the first order deviations (denoted by 1) as small com- 
pared with the background quantities (denoted by 0). The undisturbed met- 
ric 9^ is given by Q). 

Then starting from the field equations Q2.8Q , using the energy-momen- 
tum tensor ( |2.9| ) and neglecting the gas-pressure against the matter-density 
we obtain in the lowest order the well known results (Friedman equations, 
•--2-V 

dx 4 -'' 

R 2 8nG /0 ox , . e 



a 2 



(P p + Pr) + |A , (3.5) 



R 2 3c 2 V^p rr; 3 R 2 



= R = ~^ {Pv + 2 Pr) + \^ (3-6) 
and from the equations of motion (j2.3|) (see also fl3.3|)) it follows: 



p R 3 = const., p r R 4 = const.. (3.7) 



Latin coordinate indices run from 1 to 3. 
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Concerning the density fluctuations of matter and radiation given by 
p p and p T we consider only local ones in view of the fact, that the galaxy 
formation is a local phenomenon in such a sense, that for the linear size of 
the fluctuations, i.e. for their wave-length A the following relation holds: 

X = Rr<^R -> r«l. (3.8) 

Thus for the computation of the cosmic density fluctuations we can neglect 
all terms proportional to er,er 2 etc. resulting from ( |3.1| ), so that the line- 
element can be approximated by 

ds 2 = dx* 2 - R 2 (t) dx 2 + g u Jx a ) dx^dx", (3.9) 
v . ' M 

including the first order perturbations of the metric, for which we introduce 
the deDonder-gauge^ 

V> M V = o with w = g»v-\k P g aP g^ V = Vv*<K (s.io) 

After a (3+l)-decomposition the equations ( |3.10| ) read for the metric itself: 

i,9mn\n 9 9nn\m) 2 9 44 \m 9m4\4 ^ Ot 9m4 0, 

(3.11) 

\94n\n ~ 2 9nn\&) ~ 2 9 44\4 ~ & a 9 44 = U, 

where the sum convention holds with respect to Latin indices. 

Now we consider at first the density fluctuations following from ( |2.3|) in 
first order perturbation calculation. Because the fluctuation phenomenon of 
the plasma is non-relativistic (sound-waves), it is appropriate to perform a 
(3+l)-decomposition. Then we obtain from (|2.3| ) using the metric (|3.9|) and 
the decomposition ( |3.2| ) the following equations of motion and continuity 
for the plasma 

(3-12) 



l ; ^ means the covariant derivative with respect to the background metric g 
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P P (^\i + T li) +P P |4 + 3a(p p + i e !^p r ) =0 
and for the radiation 



p r (c i | 4 + rl 4 + «c i ) = -§ir 2 (i 



il 2 A 



)Pt\i + 



(3.13) 



(3.14) 



+ 



4 Pr 



V -c l ) - ±aR 2 p r g 4i 



I p r (c* H + T 4l ) + p r|4 + 4 a ( 1 - I ei^) p r = 0, (3.15) 

where terms proportional to c 2 /c 2 < 10~ 9 are neglected against comparable 
terms of the order of One; the Christoffel symbols of the first order in (|3.12|) 
to ( gig ) read: 

Hi = \R~ 2 {9u\i - 2^, 4 ), r 4i = iT 2 (o^ - \g^), (3.16) 

whereas the operator A has the form A = R~ 2 A, see Appendix A. 

The sound-wave equations for the density perturbations are obtained by 
the divergence of the equations of motion and the subsequent substitution 
of the velocity derivatives by the continuity equations. In this way we get 
with the use of the Christoffel symbols (|3.16|) and the deDonder-gauge (|3.11 ) 
from ( |3TT2D and ([TTJ) 



R 2 + Idet 1 ^ S t \i\i) - 5p| 4 |4 



2a 



4 d_ 



c- 1 7 -/ 2 A r 1 1 _ 1 

<5p|4 + 2 de9 7 ~ 0r|4 +2^14 + 2 9m\4 



(3.17) 



5 P | 4 - 1 5 r , 4 + \a (1 + I d) el ^ S r \-d(f3+f a 2 1 2 A) 5 r - 

- (« 2 + /?) (7 + 3 <? 44 ) + \ (R~ 2 g m \i ~ ^44|4| 4 ) = 
for the relative matter fluctuations 5 P and from (|3.14 ) and ( |3.15| ) 

|ir a (i- e ^)^ |4|4 _ 4 



r|4|4 



a 



+ 



(l-et^)5 r|4 + | 7 |4 + f L H 
+ I- 1 [I <5 P |4 - 5 I]4 + a (1 + § d) ei^ 5 r J + (/? + |a 2 l 7 2 A) e*'^ 5 r - 

-4/5(| 7 +^44) +l(^" 2 ^44| i | i -^44|4| 4 ) =0 (3.18) 
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for the relative radiation fluctuations 5 T with the definitions 



5, 



1,0 r 1 , ,0,0 1 „_o 1 

P P /Pp> d r = p T /p T , d = p T /p p , j = R g { 



(3.19) 



for the density contrasts, the ratio of the radiation and matter background 
densities as well as the trace of the space-like metric perturbations. 

Finally we need the field-equations for the gravitational potentials g u 
and 7. From (|2.8|) and (|2.9|) we obtain with the use of the decompositions 
(|3.2| ) and the deDonder-gauge ( p.llj ) for p, — v — 4: 



R 2 ^44|i| 4 -^44|4|4- 2 «(7|4 + § ^44 1 4 ) -2a 2 (7 + 3# 44 ) - 2 /3 7 
= ^P P (6 P + 2d8 I ) 



(3.20) 



and for the trace concerning p, v = m, n: 

i 4 - 2 « (<?44|4 + §7|4) + 6 « 2 (7 -^44) -6/3^4 + 4 

= 2 S e Pp(«s P + l^) 



D-2 1 1 



1 

7 



(3.21) 



using also the zero order equations ( p.5| ) and (|3~H|). 

The 4 coupled wave-equations fl3.17p , (|3.18| ), ( |3.20|) and ( |3.21| ) are those 
of the self-gravitating two-fluid system in consequent first order perturba- 
tion. It is very interesting that using the deDonder-gauge the wave-operator 
applied to the metric perturbations g 4i and 7 does not only appear in the 
wave-equations for the metric (|3.20|) and (|3.21|) but also in the wave-equations 
of the matter and radiation fluctuations ( p. 17] ) and ( |3.18| ), so that it can be 
eliminated there. By doing this we obtain from ( |3.17 ) and (|3.18|) simultane- 
ously substituting a 2 and (3 by ( |3.5| ) and (pij|): 

R~ 2 5 p \i\i + \de^^5 r \i\i) -5p| 4 |4-2a 5 P | 4 + \de^ 1 ^ 5 r \ A + fflW 



4 d_ 



*p|4-l*r|4 + |a(l + |d)et l 7^«J r 



+ 



+ d 



4?rG 
3c 2 



° Pp (l + 2d-f (l + rf)/ 7 2 A)-!A(l + f Z 2 A) + 



(3.22) 



+ f eiT 2 / 2 A 



+ 4 -^P p (S p + 2d5 r )+ ^p p (l + 2d)-A 



4ttG ° 
,2 



#44 = 



10 



and 



i R- 2 (1 - e WA) J r|i|i - 5 r | 4 |4 - « [(1 - e^) 5 r , 4 + |k| 



+ 



+/-■ 



■ 5 P \4 - 5t\4 + a (1 + id) e 



AttG 



3 "PI 4 

£ p (l + 2d-f (l + d)Z 2 A)-±A(l + f Z 2 A) + 

+ f eiT^Aj e§^5 r + 
+ W P p (5 P + 2d<y + f fep p (l + d)-|A 



K 2 



7 + 



+ 4 



^Pp(l + I^)-^1^44 = 0. 



(3.23) 



In earlier papers concerning galaxy formation on a general-relativistic level 
one has tried to solve the wave equations for the metric perturbation and the 
sound-wave equations for the density fluctuation simultaneously. However 
we will see now, that this is not necessary for understanding the galaxy 
formation; for this purpose the sound-wave equations (|3.22|) and (|3.23j ) for 
the matter and radiation fluctuations are totally sufficient. 



4 Newtonian limit 

The treatment of the set of the 4 first order perturbation equations ( |3.20| ) to 
(|3.23|) can be simplified significantly by transition to the Newtonian limit ne- 



glecting all terms proportional to c 2 in fl3.22| ) and ( |3.23[ ); with this approach 



the pairs of the equations ( |3.22| ), ( |3.23j ) and fl3.2C| ), (|3.21|) can be solved suc- 



cessively. From ( p.20| ) and ( |3.21| ) one sees immediately, that g u and 7 are of 



the order of c 2 and become therefore important in fl3.22j ) and ( |3.23|) only for 



fluctuations with masses M > 1O 19 M ; a detailed estimation of this value is 
given in Appendix B. Restricting ourselves in view of the galaxy masses to 
the only interesting lower mass range we can neglect in (|3.22| ) and ( |3.23 ) all 



terms proportional to g 44 , g u \ 4 and 7 (trans-Newtonian terms) |] so that the 
'Note, that the first order approximation is not identical with the Newtonian limit. 
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relevant sound-wave equations read after additionally substituting x 4 = ct: 

R~ 2 (c s 2 5 P |i|i + | dc 2 e^ A 6 t \i\i) - 6 pW t -2ac(8 p \ t + \de^ A 8 t \ t ) - 

~ I Zf [^Pl* - ! S *\t + ! « C ( 1 + I d ) el ' ? - <*r] + 

+ d[^p p (l + 2d-§ (l + d)Z 2 A) -lAc 2 (l + f Z 2 A) + (4.1) 

+ f e|,Z 2 A] e§^5 r + 

+ 4vrGp p (5 p + 2c/(5 r ) = 
and 

| c 2 Rr 2 ( 1 - eWA) 5 r | ir ^| ir «c(l- eWA) <y r|t + 
+ f [IS P \t-5 llt + ac{l + ld)el l ^S I ] - 

- [*f° Pp (l + 2d-f (l + d)/2A)-iAc 2 (l + f / 2 A) + (4.2) 

+ f e -J/ 2 A]el^5 r + 

+ i ^P p (8 p + 2d5 I ) = 0. 

It is very interesting, that these two equations for the matter and radiation 
fluctuations can be treated with the only use of the background equations 
( |3.5| ) and (ffTTp; the equations ( |3-2U| ) and ( p. 21] ) for g 44 and 7 are superfluous 
now! However these two last wave-equations for the gravitational potentials 
g 44 and 7 tell us that the matter and radiation fluctuations described by Q4.1|) 
and (|4.2j ) are the source of gravitational waves immediately correlated with 
the fluctuations (see also Appendix B and footnote 7). 

The sound-wave equations (|4. 1|) and Q4.2|) are identical with those of the 
Newtonian cosmology (GeBner & Dehnen 2000), apart from small modifica- 
tions concerning the A- and e-terms and a better consideration of the energy 
conservation in (|3.13| ) and Q3.15| ). Thus it is shown once more explicitly, that 
one reaches exactly the Newtonian limit from general relativity by using the 
deDonder-gauge (Rose & Corona- Galindo 1991). 

The integration of the equations (|4.1|) and ( |4.2| ) will be performed in the 
following within two steps. First we investigate prior to the recombination 
epoch the stable and unstable oscillation modes with the method of the dis- 
persion relation, where the unstable modes yield the mass-spectrum for the 
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galaxies. Subsequently we use these results as initial values for a numerical 
integration of the unstable modes up to the present time in order to deter- 
mine the amplitude increase of the fluctuations during the evolution of the 
Universe. The last step is required since the coefficients of the sound-wave 



equations fl4.1|) and (|4.2|) are time-dependent, which implies that the disper- 
sion relation method is only approximately valid for smaller time periods, 
see Q. 



5 The dispersion relations 
5.1 Theoretical approach 

As one can prove easily, prior to the recombination of the hydrogen at the 
redshift z R = 1380 (see ( |5.18| ) ff.) the terms proportional to A and e in the 
equations ( |3.5|) , (|3~6| ) and ( [4.1|) , (|4.2|) can be neglected in comparison with the 

p p -terms, since at the recombination time the relations hold: AttG p p / Ac 2 ~ 

1, 5 • 10 7 and 4vrG p p R 2 /c 2 ~ 3 • 10 2 (A « 1(T 56 cm" 2 ). Then we obtain from 



( |4. 1| ) and ([4.2|) with the eigenvibration ansatz 

5 V = Ae l{ujt ~ RkS \ 6 T = B e i{ - ut - R% ^ (5.1) 
with (| k\ = k = 2tt/A) 

kR = 27r(p p /M) 1/3 R = const., M = °p p A 3 = const. (5.2) 

(A wave-length and M mass of the fluctuation) the following algebraic equa- 
tions for the complex- valued amplitudes A and B: 

A [uj 2 -iuc(2a + l^) - c 2 k 2 + 4ttG ° Pp ] + 

+ B [-tooc(ade-l k21 ' - ^) - ldc 2 k 2 e-l k2L " + 8nGdp p - (5.3) 

- s^d(l + id) e"t fc % 2 + *f- p p d (l+2d+f (1 + d) k 2 l 2 ) e~l fc % 2 ] = 
and 

A[li^ + ^°P P }+ (5-4) 
+ B [to 2 - icuc {a (1 - e"t fc % 2 ) + ±) - l c 2 k 2 (1 - ) + ™fdp p + 

+ as* (i + id) e-l k % 2 - ifp p (l + 2d+f (1 + d) k 2 l 2 ) e"l fc % 2 ] = 0. 
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However, since the coefficients in ( |5.3| ) and ( |5.4j ) (see also ( |4.1| ) and ( |4.2| )) are 
time-dependent due to the cosmic expansion, the frequencies u become time- 
dependent. Only if the coefficients in and ( |5.4| ) change slowly compared 
to the age of the Universe, then also u changes slowly with time as a param- 
eter, and the procedure with a nearly constant uj is a useful approximation 
method for not too long times. In other words: Deriving (|5.3|) and ( |5.4|) the 
time derivative of 5 PjI ~ e lujt reads strictly i(u + ujt)e tuJt , but the second 
term within the bracket has been neglected under the assumption 

\u\t < |a>|. (5.5) 

This condition for the applicability of the dispersion relation method will be 
proved retrospectively. 

The two linear homogeneous equations (|5.3| ) and (|5.4| ) for the amplitudes 
A and B only have a non-trivial solution, if the determinant of its coefficient 
matrix vanishes. This condition results in a fourth-order equation for the 
frequency u(k) (dispersion relation) 

u A + iLU 3 A-LU 2 B-iujC + D = (5.6) 

with the (time-dependent) coefficients: 

A = -3ac(l-le-l k2l t)-±(l + ld), (5.7) 

B = |c 2 A; 2 (1 -e-i fc % 2 ) + 2 ac [ac (l - e"§ fc % 2 ) +JL(l + |d)] - 

-4vrGp p [l + |d-|(l + 2rf+f (l + d)A: 2 / 2 )e-t fc % 2 ] - 
-^(l + f d)e-t fc % 2 , (5.8) 



C = ac 



+ 



- | ( 1 + 2 d + f (1 + d) k 2 I 2 ) e~i fc2 ' 2 ] - 2^ (1 + | d) e~l k ' 
+ ± [4vrGp p (l + 4rf+f d 2 )-ldc 2 k 2 }, (5.9) 
D = c 2 k 2 [\c 2 k 2 {\- e"l fc2 ' 2 ) - eg (1 + | d) e~l k % 2 ] - 

-^G° Pp [ c 2 k 2{l-{l + \d) e-§ fe % 2 ) - 3 ( 1 + | d) 2 eg e-f fe % 2 + 
+ 47rGp p (l + |rf)(l + 2d+f (l + d)k 2 l 2 )e-l k21 *]. (5.10) 
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In these expressions we have neglected again c 2 k 2 - against comparable c 2 k 2 - 
terms. With the solutions u>(k) of (|5.6|) the amplitude ratio r(k) = A/B is 
determined by ( (5.3| ) or ( |5.4| ), which is in general complex- valued including 
the phase relation between matter and radiation fluctuations according to 

7" — | 7" | 

Since the frequency u as solution of ( |5.6| ) is complex- valued, we set: 

u(k)=a(k)+ib(k), a(k), b(k) e K ; (5.11) 

here a(fc) represents a real-valued frequency connected with the period T = 
27r/| a | and r = denotes a real- valued relaxation time of the fluctuations, 



so that in view of ( |5.1| ) b < describes increasing and b > decreasing 
amplitudes, i.e. r < means growth times and r > decay times. Thus 
unstable modes are uniquely characterized by 

a = 0, 6<0, (5.12) 

which give rise to a kind of spinodal decomposition of the homogeneous 
mixture of plasma and radiation. 

Insertion of (|5.11|) into (|5.6 ) and separation of real and imaginary parts 



result in the two coupled relations for a(k) and b{k): 

a 4 - a 2 {6b 2 + 3 Ab + B) + 6 4 + Ab 3 + B b 2 + C b + D = 0, (5.13) 

a [a 2 (Ab + A) - Ab 3 - 3 Ab 2 - 2 B b - C] =0. (5.14) 

Because we will restrict ourselves in view of the galaxy formation to the un- 
stable modes defined by ( |5.12j ) equation (|5.14|) is fulfilled identically, whereas 



equation (|5.13|) is reduced to: 

b 4 + Ab 3 + Bb 2 + Cb + D = 0, (5.15) 

the negative solutions of which yield the growth times as a function of k 
or M according to k = 2ir (pp/M) 1 ^ 3 (see (p.2|)). A detailed discussion of 
all solutions of (|5.13|) and (|5.14|) is given in our previous paper (Gefiner & 



Dehnen 2000). For the ratio r = A/B we deduce from (f)]4|) in the case of 
a = 0, see (|5l2D : 

6 2_ 6c ( a(1 _ e -|fc% 2 ) + ] L) + i ( a k 2( 1 _ e -%k*i2 ) _ ^G d p p ~ 
- eg (1 + Id) e-i fc % 2 + ±f° Pp (1 + 2d+ f (1 + d) k 2 l 2 ) e~l k21 
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Evidently r is real-valued, i.e. the nonoscillating fluctuations of matter and 
radiation are in phase (r > 0, acoustic mode) or in phase opposition (r < 0, 
optical mode). 

For the numerical evaluation of ( |5.15|) and ( 5.16|) the knowledge of the 
coefficients A,B,C and D prior to the hydrogen recombination is required. 
To this end we reduce the quantities in (|5.7p to ( |5.10| ) with the help of the 
"redshift"^ z = (R H /R) to the today's background values denoted by the 
index H: 



(5.17) 

2 — 2 rl — rl °_°3 

c s ~~ c sh z i a — a H z, p p — p pH z . 

In this context we take for the present background matter density and the 
present background radiation temperature the values 

p pH = 2-10- 31 g/cm 3 , f H = 2,7K, (5.18) 

where the matter density corresponds to the baryonic one determined by the 
abundance of the light chemical elements generated in the cosmic nucleo- 
synthesis. Accordingly, the hydrogen recombination has taken place with 
regard to the Saha-equation at the redshift z = z R = 1380 and the today's 
background quantities take the values (see (|2.5| ) and (|2.6|) ): 



= P rH /p pH = aT 4 H /(p pH c 2 ) = 2,24-10- 3 



'pH " -*■ H / vrpH 

. -, ... , , V HJ 



Z 7H = m/(a Th p pH ) = 1,26- 10 cm, 
c s 2 H = 2 k B T H /m = 4, 46 • 10 8 cm 2 /sec 



(5.19) 



k H = 2 7 r(p pH /M) 1/3 = 2,92-10- 21 iV- 1 / 3 



cm 



where N is the number of solar masses (M = NM & ) of the matter fluctua- 
tions. 

With these data we have to solve the fourth-order equation ( |5.15| ) for 
b, which in general yields either 4 real-valued solutions, or 2 real-valued 
solutions and 2 complex ones, or 4 complex-valued solutions; however, only 
the real-valued solutions are relevant for our purpose. Thus we seek, with the 
help of a standard computer algebra program, these solutions as functions of 



3 Strictly speaking it holds z = (redshift + 1) 
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N for several values of z or s = z R /z around and prior to the recombination 
of the hydrogen [z ~ z R = 1380, s ~ 1). Subsequently we calculate the real- 
valued amplitude ratio Q5.16| ) also as function of N. The results are shown 
in the Figs. 1-4. Finally we prove the condition (|T|), which reads for b 
using (|6.6|) and setting z = z R j s 



db/ds 



< \s/ (s + d H z R )(s-2d H z R ) + 2(d H z R ) 3 / 2 y/s + d n 



(5.20) 



for s < 1. The left hand side of Q5.20| ) can be rewritten with the help of 
(1515D as 



db/ds 



dA^ + dB b2+ dC b+ dD 
ds ds ds ds 



Ab 3 + 2Bb 2 + 3Cb + AD 



(5.21) 



and subsequently calculated for fixed N at every point (time/redshift) s by 
the use of (O) to <KW) and flSTTTD. 



5.2 Results 

Fig. 1 shows qualitative courses of the solutions for b along with the ratio 
r in the case of a = (nonoscillating modes) shortly before the hydrogen 
recombination. The mode b\ is a purely decaying optical one (6 > 0, r < 0) 
and has decay times r = b^ 1 of the order of 2,3 ■ 10 3 y > r > 1,8 • 10 3 y, 
which is shown in Fig. 2, so that this mode is damped out very rapidly at 
the recombination time of 4, 8 • 10 5 y (see ( |6.6| ) ff.). Of special interest is, 
however, the mode b<i : first, a mode transition from an optical to an acoustic 
mode lying in the damped range (b > 0) occurs at 2, 01 • 10 5 M , so that 
at this point the radiation fluctuation mode vanishes (|r| — > oo). Secondly 
this mode has an unstable acoustic region {b < 0, r > 0) for 2, 3 • 10 6 M Q < 
M < 5,8 • 10 12 M Q , corresponding to scales between 219 ly and 2,9 • 10 4 ly 
just before the recombination time. The ratio r in this unstable range is 
very large, i.e. of the order of 10 6 up to 10 8 (see also Fig. 6), wherefore the 
radiation fluctuations are very small compared to the matter fluctuations, so 
that it is questionable, whether an observable influence on the background 
radiation data exists.[] 

7 The radiation density contrast is correlated with a temperature fluctuation of the 
background radiation according to <5 r = 4AT/T. A second contribution to the tempera- 
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This instability region coincides exactly with that of the observed masses 
of galaxies including globular clusters. Fig. 3 shows the behaviour of b and 
t = 1/6 in detail for the mode b 2 at s = 1; the growth time is of the order of 
10 9 y, so that the unstable growing mode needs an additional strong increase 
during the later evolution of the Universe for galaxy formation (see also Bon- 
nor 1956), which requires a positive cosmological constant (quintessence) and 
positive curvature as discussed in the next section. The whole instability val- 
ley as function of "time" s is illustrated in Fig. 4. The lower mass limit is 
the Jeans mass, where gravity exceeds the plasma pressure; the upper mass 
limit has its origin in the increasing effective radiation pressure on the plasma 
with increasing size of the fluctuations. The lower and the upper mass limit 
coincide at z ~ 7600, so that the instability valley ends with a mass value 
of 2, 3 • 10 6 M at the time of 1, 7 • 10 4 y (s = 0, 18; T = 20520 K). This fits 
approximately to that time/redshift (z ~ [m? G / T( kT H p m ] l / A a T ^ 2 ), where 
the free path length of the photons (upper mass limit) and the Jeans length 
(lower mass limit) become identical, and which consequently represents the 
earliest time for galaxy formation! All masses within the valley show in- 
creasing instability with time and also remain unstable after the hydrogen 
recombination. Their further time evolution will be investigated in the next 
section. 

Finally Fig. 5 shows the applicability condition (|5.20|) of the dispersion 
relation for the modes b\ and b 2 and typical masses versus s. In the case of 
the unstable growing mode b 2 the condition is fulfilled more or less, so that 
referring to this our results presented above are decidedly useful. 

6 Time evolution of the unstable fluctuation 
modes 

6.1 Theoretical approach 

Since the dispersion relation method is valid only locally compared to the 
timescale of the cosmic expansion, we investigate the exact time evolution of 

ture fluctuation results in consequence of redshift from the variation of the gravitational 
potential according to AT/T = g A4 /2, see equ. ( |B.3| ). 
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the density contrasts with the ansatz (instead of (|5.1|)) 



f(t)e 



—iRkx 



5 r = g(t) e 



iRkx 



(Rk = const .) 



(6.1) 



and obtain from (|4.1|) and (|4.2| ) the exact differential equations for the am- 
plitudes f(t) and g(t) ('= J^): 

/+ [2ac+ |d^] f + d(ae~l l ' k2 - g + c 2 k 2 f + 

e I fc 3 +^(l + |(i)-fJ p (l + 2d+f(l + (i)i J l 7 J ) + 



, _4,2r.2 
+ 0te 9 i 



-4nGp p (f + 2dg) = 
and 

<? + [ac(l-e-l k2l t) + f]g-lf- f + l c 2 k 2 g- 



9 ~ 



(6.2) 



_i 1,2,2 
g 9 7 



ic 2 fc 2 + ^(l + |ci) 



4ttG ° 



p p (l + 2rf+f (l + rf)A; 2 / 2 ) + 



+ |Ac 2 (l-f P/ 2 ) + f ikH- 



g- 



WnG 



p(f + 2dg) = 0. 



(6.3) 



Additionally the differential equations ( |3.13| ) and fl3.15| ) can be integrated 



immediately with the use of (|6.1|) ; one gets for the flow- velocities in the 
Newtonian limit 



(f + adge l 1 ^ 2 ) 



k % e 



iRkx 



iRk 2 



c l = l(g-age ^ k ) 



i „—i Rk x 



k l e 



iRk 2 



(6.4) 



from which the real or imaginary part is to be taken correspondingly to fl6.1| ). 

Because the coefficients in ( |6.2|) and flfOl ) are known as functions of z in 
view of (|5.17|) we substitute the time t through z as an independent variable 
according to (c.f. (|3.5| )): 



dz 

z 



-acdt 



^p pH ^ 3 (l + ^) + |Ac 2 -^^ 2 



1/2 



dt. 



(6.5) 
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Simultaneously we have taken into account the A- and e-terms in 
and (|6.5|) and also in a given by ( |3.5|) , in order to integrate over a timescale 
starting prior to the recombination epoch (z > z K ) up to the present time 
(z = 1). In this way the influence of the cosmo logical constant A and the 
curvature sign e on the time evolution of the amplitudes of the (unstable) 
fluctuations will be investigated. 

The differential equation (|6.5|) for z(t) cannot be integrated exactly, so 
that the integration of (|6.2j ) and (|6.3|) with respect to z is even necessary. 
Only for the early stages of the Universe, where the A- and e-terms on the 
right hand side of ( |6.5| ) can be neglected, i.e. for z > z R , the integral can be 
given explicitly as follows: 

(=i(¥U ^ 3/2 [ VTTd^ (l-2d ll z) + 2(d ll z)^}. (6.6) 



According to this solution (see also (|5.18 ), ( 5.19|) ) the recombination of the 
hydrogen (z R = 1380) has happened at t — 4, 79 • 10 5 y. 



In the case, that the eqs. Q6.2|) and ( |6.3|) are integrated over a timescale 
across the hydrogen recombination, we have to consider the decoupling of 
matter and radiation due to the hydrogen recombination. This means in 
detail that for z < z K holds (instead of ( |5.17| )) 



assuming an adiabatic cooling and an adiabatic sound-velocity for the de- 
coupled neutral hydrogen. In this context we have determined the detailed 
behaviour of Z 7 (and Z e ) around z = z n with the help of the actual ionisation 
degree according to the Saha-equation. 



For the numerical integration of (6.2) and (|6.3|) after consideration of (|6.5| 



and (|6.7| ) we take as initial values those of the (unstable) nonoscillating eigen- 
modes given by the dispersion relation prior to the hydrogen recombination. 
For this purpose we identify (|5.1|) with ( |6.1|) setting: 



f = Ae iujt , g = Be iuJt (6.8) 
and get with Aj B = r = \r\e l< ^ 

f =\r\e i{Lut+(t '\ g = e iujt (6.9) 
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up to an arbitrary common constant factor. The complex-valued ansatz 



(6.8) for solving the set of the linear homogeneous differential equations ( |6.2|) 
and ( |6.3j ) for the real-valued density contrasts means, that the real and the 
imaginary part of ( |6.9| ) constitute the two independent real-valued solutions 
of the system. 

Considering at first the imaginary part we obtain after substitution of uj 
according to (|5.11|) 



/ = \r | e~ bt sin(at + 0), g = e~ bt sinat, 

/ = -\ r \be~ bt sin(at + 0) + \r\ae~ bt cos(at + 0), (6.10) 

g = —b e~ bt sin at + ae~ bt cos a t. 

Now for the nonoscillating modes with a = the amplitude ratio r is real- 
valued (see ( |5.16| )), i.e. = 0,7r. Herewith it follows from ( |6.10| ) 

/ = 0, (7 = 0, / = 0, g = (6.11) 

for every initial time t, so that the imaginary part leads only to the trivial 
solution. 

In the case of the real part it follows from ( |6.9p 

/ = \r | e~ bt cos(at + 0), g = e~ bt cos at, 

f = -\ r \be~ bt cos(at + 0) - \r\ae~ bt sin(at + 0), (6.12) 

g = —be~ bt cos at — ae~ bt sin at. 

Restricting ourselves again to the nonoscillating modes with a = we find 
for the acoustic modes (0 = 0) 

f=\r\, g = l, f = -\r\b, g = -b (6.13) 

and for the optical modes (<ft = 7r) 

f=\r\, g = -l, f = -\r\b, g = b (6.14) 

as initial condition for every initial time t unique up to a common constant 
factor. The values of the 'frequency' b and the amplitude ratio | r | for the 
initial time t connected with the initial "redshift" z according to ( |6.6|) must 
be taken from the dispersion relation (see ( |5.15| ), ( |5.16| ) and Fig. 1) after 
choosing the mass M of the fluctuation. Then the integration of ( |6.2j) and 
(6.3) can be performed numerically starting before the hydrogen recombina- 
tion at z > z R (= 1380), where for technical reasons it is appropriate to use 
instead of z the new independent variable s = z R J z. 



21 



6.2 Results and Conclusions 



The foregoing description of the numerical integration of the differential equa- 
tions ( |6.2|) and (|6.3| ) for the fluctuation amplitudes has been performed for 
several cosmological models. In all cases we have chosen t — 1, 25-10 5 y = s = 

0. 5 as starting time. The unstable modes are all acoustic ones (see Fig. 1), 

1. e. the initial condition fl6.13| ) holds. As masses we choose N = 4 ■ 10 6 and 
N = 7, 9 • 10 9 , which represent the boundary masses of the instability valley 
for the initial time; from the dispersion relation we then get the values for 
b and |r|, see Fig. 6. Fig. 7 shows the behaviour of the matter fluctuation 
amplitudes f(s) up to the present time (s = 1380) for several situations: 

• case a) for A = 2 • 10~ 56 cm" 2 , e = 0; 

• case b) for A = 0, e = 0; 

• case c) for A = 3, 279 • 10~ 56 cm" 2 , e = +1, R n = 3, 388 • 10 28 cm, 

H = 93,4km/secMpc; 

• case d) for A = 2, 206 • lO" 56 cm" 2 , e = +1, R H = 3, 306 • 10 28 cm, 

H = 75km/sec Mpc; 

for the matter density holds (|5.18|) in all cases. Case c) is identical with a 
cosmological model proposed by W. Priester et al. (1995). Case d) represents 
the critical one, where i? = 0,-R = 0ats = 282 (redshift 3,9). The last two 
cases c) and d) show clearly, that a cosmological constant (or quintessence) 
together with e = +1 promotes the galaxy formation very strongly by an 
increase of the unstable matter density contrast up to 4 orders of magnitude, 
whereas in the two first cases a) and b) the increase is too small for galaxy 
formation. In this context it may be of interest, that the strong increase of 
the matter density contrast (region with a ~ 0, see below) occurs at red- 
shifts of the order of those of the quasars. On the other hand, a cosmological 
constant or (nonbaryonic) cold dark matter alone, that means an increase 
of a (c.f. (|3.5| )), reduces the tendency for galaxy formation drastically, see 
case a ). The physical reason for this feature lies in the damping terms (a f- 
and ag-terms) of the differential equations (|6\2] ) and (|6.3|). Only if these 
terms go down for a longer time, and this happens according to ([O]) only 
by a compensation through the e-term with e = +1, then the increase of the 
density contrast of matter approaches to an exponential one with respect to 
time. These results favour a cosmological model with a positive cosmological 
constant and a positive curvature (Lemaitre-Universe), as also proposed by 
W. Priester et al. (1995) - however for other reasons - and are in agreement 
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with the analysis of the angular power spectrum of the background radia- 
tion (McGaugh 2000; Overduin & Priester 2001). Of course, this statement 
contradicts the conjecture of a flat Universe concluded from the inflationary 
cosmological model (c.f. the final remark). On the other hand, a positive 
cosmological constant is also in accordance with the observation of an accel- 
eration of the cosmic expansion at the present time (Riess et al. 1998). 

Additionally, the behaviour of the density contrast g(s) of the radiation is 
given in Fig. 8 for the same cases as in Fig. 7; after the decoupling from the 
matter due to the hydrogen recombination, the radiation fluctuations decay 
finally very rapidly as expected for the considered sizes. However, only in 
the region of a very strong increase of the matter density contrast also the 
radiation density contrast increases slightly in consequence of gravitational 
coupling. But this property is no longer that of an eigenmode; obviously the 
unstable matter-radiation system leaves after some time the initial eigenmode 
in consequence of the cosmic expansion. 

In order to confirm the foregoing discussion of the influence of nonbaryonic 
cold dark matter coupling only gravitationally the calculations presented 
above are repeated, after in the cosmological background equations (|3.5|) and 



(|3.6f) the matter density p p has been enlarged by addition of several amounts 

of cold dark matter as follows: p p — > p p (1 + 5) with 5 = 0,1; 0,25; 0,5; 1; 2,5; 
5; 7,5; 10.0 The results for the matter and radiation fluctuation amplitudes 
are illustrated in Fig. 9 and Fig. 10 for the critical case (with R = R = 0), 
for which we have obtained in Fig. 7 (case d) the strongest increase of the 
matter fluctuation amplitudes; the initial time is again s = 0, 5 and as mass 
we choose M = 1, 8 • 10 8 M , which represents the mean value of the unstable 
masses at the initial time. The comparison with the case 5 = shows, that 
even a small amount of nonbaryonic cold dark matter reduces the increase of 
the matter density fluctuations drastically, so that in view of the existence of 
galaxies no essential amount of nonbaryonic cold dark matter may be present 
in the Universe, which coincides also with other new results, for instance with 
the analysis of the angular power spectrum of the background radiation (for 
a review see Overduin & Priester 2001). 



8 It is assumed, that the dark matter does not participate in the fluctuations. Because 
of the gravitational coupling this assumption is only approximately true. 
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7 Final Remark 



The results of our calculations compared with the observations of galax- 
ies in the mass range of 2 • 1O 6 M < M < 6 • 1O 12 M favour a cos- 
mological Robertson- Walker model with a positive cosmological constant 
(A ~ 2, 2- 1CT 56 cm~ 2 ) and a positive curvature (Lemaitre-Universe), however 
without an essential amount of nonbaryonic cold dark matter. This is very 
satisfactory for two reasons: First of all, the age of the Universe will be of 
the order of approximately 30 ■ 10 9 y (c.f. Priester et al. 1995), so that no 
contradiction exists with respect to the age of the oldest stars. Secondly, a 
positive curvature connected with a finite (spherical) space is the only case, 
where the properties of the Universe as a whole, as e.g. homogeneity and 
isotropy can be proved empirically up to a certain extent. In all other cases 
this is totally impossible, wherefore the cosmological models with vanishing 
and negative (hyperbolic) curvature, i.e. infinite space, lie beyond the scope 
of natural sciences. Consequently, the interpretation of the results of the 
inflationary model as necessity for flatness of the Universe as a whole is very 
problematic and not convincing (see also Ellis 1988; Hiibner & Ehlers 1991). 
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Appendix A 

Derivation of the radiation pressure in an op- 
tically thin medium 

We consider photons confined by a plasma. The usual radiation pressure p r 
on the plasma is then proportional to the radiation energy density u = p v c 2 : 

p x = \u. (A.l) 

However, this is only correct if the plasma is optically thick, so that u remains 
nearly constant with time. This does not apply in an optically thin plasma 
because of the stronger diffusion of the photons and their energy. 

Thus we have to derive in a first step the relativistic diffusion equation 
for the radiation energy. For this we start from the energy-momentum tensor 
for radiation [c il &' = 1) 

T/ = uc^-p r (<V-^) (A.2) 
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and obtain with ( |A.1| ) and T^ U ]W = and subsequent multiplication with c M : 

{uc,\ w g^ = \u\^. (A.3) 

With regard to the decompositions ( |3.2|) up to ( |3.4|) equation (|A.3|) is ful- 
filled in the lowest order; but in the first order we obtain after a (3+1)- 
decomposition 

f (uc m )| n £™ + U|4 + 4au = 0, (A.4) 

where all terms proportional to c~ 2 are neglected simultaneously (Newtonian 
limit). With the assumption of Fick's law for the diffusion flow 

uc m = -qIjUIui (A. 5) 



(A 7 free path length of the photons) and assuming a rapid diffusion compared 
with the cosmic expansion we get from (|A.4 ) the following diffusion equation 
for the radiation energy: 

R~ 2 u\ m \ m - u\ t = 0. (A. 6) 

Now, the effective radiation pressure p 7 ^ e on the plasma is determined by 
the photon energy density remaining after a mean collision time r = Z 7 /c of 
the photons with the plasma electrons. Thus we obtain 

P^M = ±u(t + r), r<t (A.7) 
with the Taylor expansion 

P^e(t) = lEZoh^ n) (t)r n . (A.8) 
On the other hand, from the diffusion equation (|A.6|) it follows by iteration 



and under the assumption, that R(t) ~ R(t + r): 

£(»)(*) = (| c / 7 i?- 2 A)"M(t). (A.9) 
Insertion of ( |A.9|) into ( |A.8| ) gives with r = l 7 /c: 

Pr-(*) = lE^oA(l# A ) B «(*) = 

(A.10) 

= | el^ u(t) = el ^ p r , A = R- 2 A. 



26 



The same consideration is valid for usual nonrelativistic gas particles, where 
the diffusion equation for the particle density is identical with (|A.6|) after 
substituting the factor 9/4 by 3. The operator e^ l ~<— in ( |A.10| ) takes into 
account that the pressure is determined by the particles (photons) remaining 
after a mean collision time and that it takes effect only on scales larger than 
the free path length / of the particles under consideration. 



Appendix B 

Estimation of the order of magnitude of the 
trans-Newtonian g u - and 7-terms 

The differential equation ( |3.2(J| ) for the gravitational potential g M reads in 
its leading terms with respect to c -2 , i.e. in the near zone (Newtonian limit): 

R- 2 g ullll = 8 -¥°p p (5 p + 2d6 T ). (B.l) 

A similar Poisson-equation is valid for 7, wherefore we restrict ourselves to 
the discussion of the g 44 -terms. With the ansatz analogous to (|5.1| ), (|6TT|): 

g M = C e *M-***) , g u = h{t) e~ iR%s (B.2) 

we obtain from ( |B.1| ) immediately the gravitational potential: 

g M = -^p p (S p + 2d5 r ). (B.3) 

A comparison with the equations (|3.22|) , ( |3.23| ) shows now, that the g 44 -terms 
can be neglected so long as 

m K = B hn /3 M 2 J 3 z « 1 (B.4) 

is valid. For the hydrogen recombination epoch (z R = 1380) one finds with 
the data ( |5 . 1 8| ) 6, 8 • 10 19 for neglecting the trans- Newtonian terms in 

( |3.22|) , ( p.23| ). For decreasing z- values the critical N- value even increases. 
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Figure Captions 



Fig. 1: Qualitative courses of the frequency b (a) and the amplitude ratio 
r (b) for the nonoscillating modes (a = 0) versus \gN. There exist two 
main modes 61 and 62, where the sign characterizes the optical (— ) and the 
acoustic (+) ones. The numbers refer to the recombination epoch (s = 1). 
Fig. 2: Decay time r of the optical mode 61 versus lg N at the recombination 
epoch (s = 1). 

Fig. 3: Behaviour of the frequency b and the growth time r = b^ 1 for the 
acoustic mode 62 at s = 1. Maximum unstable region: 2, 34 • 1O 6 M < M < 
5,75-lO 12 M . 

Fig. 4: Instability valley: The negative b- values of mode 62 versus lg N and 

time s. The valley ends at s = 0, 18 and N = 2, 3 • 10 6 . 

Fig. 5: Dispersion relation condition ( |5.20|) of the modes b\ and 62 for the 

typical masses AT = 3, 16 • 10 7 (dotted line) and iV = 10 9 (solid line) in 

dependence of s. The condition is fulfilled if the numerical value of dc is 

smaller than one, that means in the case of the mode 62 for the unstable 

region. 

Fig. 6: The mode 62 and the associated values of the amplitude ratio r for 
s = 0, 5 in dependence of IgN. 

Fig. 7: The behaviour of the matter fluctuation amplitudes f(s) normalized 
to its initial value /, as function of s for N = 4 ■ 10 6 (dotted line) and iV = 
7, 9 • 10 9 (solid line), which represent the boundary masses of the instability 
valley at the starting point s = 0, 5 (see Fig. 6), in the case of several 
cosmological models: a) A = 1 • 1CT 56 cm' 2 , e = 0; b) A = 0, e = 0; c) 
A = 3, 279 • 1(T 56 cm" 2 , e = +1, d) A = 2, 206 • 10" 56 cm" 2 , e = +1. 
Fig. 8: The behaviour of radiation fluctuation amplitudes g(s) normalized 
to its initial value gi as function of s for the same cases as in Fig. 7. 
Fig. 9: The behaviour of the matter fluctuation amplitudes f(s) normal- 
ized to its initial value ft in the case of the existence of several amounts of 
nonbaryonic cold dark matter; the numbers denote the ratio 5 of the non- 
baryonic and the baryonic matter density. Starting point s = 0, 5 and mass 
N = 1,8- 10 8 . 

Fig. 10: The behaviour of the radiation fluctuation amplitudes g(s) nor- 
malized to its initial value gi for the same cases as in Fig. 9. 
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